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Coherent optical bichromatic forces have been shown to be effective tools for rapidly slowing and 
cooling simple atomic systems. While previous estimates suggest that these forces may also be ef¬ 
fective for rapidly decelerating molecules or complex atoms, a quantitative treatment for multilevel 
systems has been lacking. We describe detailed numerical modeling of bichromatic forces by direct 
numerical solution for the time-dependent density matrix in the rotating-wave approximation. We 
describe both the general phenomenology of an arbitrary few-level system and the specific require¬ 
ments for slowing and cooling on a many-level transition in calcium monofluoride (CaF), one of the 
molecules of greatest current experimental interest. We show that it should be possible to decelerate 
a cryogenic buffer-gas-cooled beam of CaF nearly to rest without a repumping laser and within a 
longitudinal distance of about 1 cm. We also compare a full 16-level simulation for the CaF B -H- X 
system with a simplified numerical model and with a semiquantitative estimate based on 2-level 
systems. The simplified model performs nearly as well as the complete version, whereas the 2-level 
model is useful for making order-of-magnitude estimates, but nothing more. 

PACS numbers: 37.10.Mn, 37.20.-l-j, 02.60.Cb 


I. INTRODUCTION 

During the past few years, methods for direct laser 
slowing and cooling of small molecules have progressed 
remarkably, advancing from a demonstration of modest 
forces in 2009 all the way to a full realization of magneto¬ 
optical trapping at ultracold temperatures M- Never¬ 
theless, all-optical slowing and cooling is presently ap¬ 
plicable only to a small number of molecules, and the 
present schemes require complex optical configurations 
with multiple repumping beams to recycle atoms that 
are lost by unwanted radiative decay into distant “dark” 
states. In 2011 we proposed that the coherent optical 
bichromatic force could improve greatly on this situation, 
by providing not only a much stronger force, but one that 
allows much greater momentum transfer prior to radia¬ 
tive loss into non-cycling dark states 0. Our 2011 paper 
relied largely on semi-quantitative arguments based on 
a two-level model, admittedly not fully adequate to de¬ 
scribe the multilevel manifolds inevitably encountered in 
real molecules. Here we verify and greatly extend this 
treatment by describing detailed numerical modeling of 
the bichromatic force (BCE) in realistic multilevel sys¬ 
tems. 

The bichromatic force in two-level systems has success¬ 
fully been demonstrated in several atoms and is amply 
described in prior publications [M||. In brief, the BCE 
is the coherent force produced by a balanced counter- 
propagating pair of two-color cw laser beams. Each two- 
color beam can be regarded as a beat note train produced 
by beams with frequencies to±6, for a carrier frequency w 
centered close to the two-level resonant frequency. If the 
laser power is adjusted so that each beat is approximately 
a TT-pulse, and if the counterpropagating beat notes are 
sequenced properly, a large net force is produced by alter¬ 
nating cycles of coherent excitation and stimulated emis¬ 
sion. Larger bichromatic detunings S produce a faster 


beat note period and thus a larger force, at the cost of 
a laser power requirement that scales as S^. In practice 
detunings of 100-400 times the radiative decay lifetime 
are workable, producing average forces larger than the 
radiative force by a similar multiple. The velocity range 
of the force is also superior to the radiative force, extend¬ 
ing over a range of S/k, where k is the wave vector for 
the transition of interest. Recent developments in diode 
laser technology and rf electronics have made it relatively 
easy to produce the multi-frequency beams required for 
the BCE. 

In a molecule the maximum attainable velocity reduc¬ 
tion will usually be limited not by the velocity range 
of the force, but instead by unwanted radiative decay 
to non-cycling dark states, typically vibrational or ro¬ 
tational levels that are addressed neither by the BCE 
laser nor by added repumping lasers 0. Thus there is a 
dual motivation for using stimulated optical forces when 
working with molecules: to enhance the average force 
and more importantly, to increase the number of opti¬ 
cal cycles that occur before the molecule is lost into an 
inaccessible vibrational or rotational level. 

In a multi-level system like the one depicted in Eig. 
[2 it is impossible to achieve a tt pulse for every compo¬ 
nent of a transition, even for a stationary molecule with 
an optimal laser pulse. Thus the optimal laser power 
and detuning involve unavoidable compromises, and the 
force is an average over the effects of several concurrently 
cyling transitions, some of them coupled coherently. In 
our 2011 paper we treated this problem phenomenolog¬ 
ically for the CaF molecule an approach repeated 
more recently for MgF by Dai and coworkers [T^. By 
contrast, in this paper we describe detailed numerical 
calculations that fully incorporate multilevel coherence, 
optical pumping, and the effects of repumping lasers that 
can optionally be used to repopulate the cycling system 
from dark states. We begin with some general predictions 
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FIG. 1. A schematic diagram of the general multilevel system 
simulated here. Ground (g) and excited (e) levels with multi¬ 
plicities Mg and Me are coupled by a bichromatic laser field 
with frequencies cu^g ± 5, and excited levels are coupled to a 
distant (d) level by a laser field with frequency There is 
no direct coupling or decay between the d and g states. Levels 
are labeled by their characteristic frequencies. 


and observations for three-level and many-level systems, 
then focus specifically on the CaF molecule, for which 
experimental tests are currently underway in our labora¬ 
tory. 

In Section|n]we describe the theoretical approach used 
to numerically solve for the time-dependent density ma¬ 
trix of a multilevel system in a pair of bichromatic laser 
fields in the rotating-wave approximation. Direct numer¬ 
ical integration has proven to be surprisingly stable and 
rapidly convergent, and we have successfully treated sys¬ 
tems with up to 16 levels. We start by treating three-level 
lambda-type systems in Section HIIl then introduce the ef¬ 
fects of dark-state repumping from distant dark states in 
Section |lVl In Section |V] we consider the additional com¬ 
plication of intra-system dark states, which can arise in 
many-level systems with high ground-state multiplicity. 
Finally, in Section IVTl we consider in detail the feasibility 
of bichromatic slowing for the CaF molecule, a species of 
considerable current interest for laser slowing and cooling 
Q. We treat the CaF problem using both full 16-level 
simulations for the B ^ X (0 — 0) Pii(1.5)/^(5i2(0.5) 
branch and a more approximate approach, in which the 
system is simplified by dividing it into several subsystems 
that are assumed to cycle nearly independently. We show 
that the subsystem-based modeling yields results very 
similar to the exact treatment in most cases, supporting 
the use of similar simplifications for treating the bichro¬ 
matic force in other molecular systems. More generally, 
the results verify the predictions of the semi-quantitative 
reasoning in our 2011 paper indicating that BCF slowing 
of molecular beams shows great promise. 


II. THEORY AND METHODS 


In our modeling of multilevel systems, we assume illu¬ 
mination by a pair of counterpropagating two-frequency 
laser beams that provide equal electric field magnitudes 
in all four components. In each beam the two frequen¬ 
cies uj + S and uj — 6 produce a train of beat notes, as 
described above. The total on-axis electric field for this 
arrangement is approximately [l^ 

EBCF{z,t) = 4 EoRe[{e{cos{kz) cos{6t) cos(x/2) 

-I- i sin(fc 2 :) sm{St) sin(x/2)}e“*‘^*], (1) 

valid so long as the beat length c/d is much longer than 
the length scale over which the system is evaluated. In 
this approximation the phase difference x between the 
electric fields of the counterpropagating beat notes re¬ 
mains approximately constant. Eq is the amplitude of 
the electric field of each of the four beams and e is the 
polarization, assumed to be identical for all four compo¬ 
nents. 

When this field is incident on an atom or molecule with 
transitions \j) O |i), each such transition is associated 
with a complex Rabi frequency defined by 




(f| d- E{z,t) \j) 
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( 2 ) 


where d is the electric dipole operator. It is important 
to note that systems may have closely spaced transitions 
with differing dipole transition elements, so the same op¬ 
tical field will simultaneously drive multiple transitions 
at different Rabi frequencies. We separate each of these 
Rabi frequencies into two terms with amplitudes ilij 
and il*j that describe co-rotating and counter-rotating 
components in the rotating wave approximation; 

= + (3) 

It is also often convenient to refer to the “Rabi fre¬ 
quency amplitude” of a transition, which we define as 

(4) 

As shown in Fig. [1] our simulations assume a system 
with a lower-state manifold of Mg levels and an upper- 
state manifold of Me levels, for which each of the lower- 
state levels may be coupled to any of the upper-state 
levels by the BCF optical field, depending on selection 
rules and line strengths. A single distant dark state is 
at times included as well, where “distant” in this setting 
means that couplings of this state to any of the others 
by the primary BCF optical field are negligible. These 
systems will be referred to in this work as Mg + Me, with 
(or without) a distant dark state. Sums over a ground 
state index (p, q, r) should be understood to run from 
1 to Mg, and sums over an excited state index (i, j, k) 
from 1 to Me- 
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The general BCF Hamiltonian includes the zero-field 
energies iJg and the electric dipole coupling due to the 
applied optical field. This Hamiltonian is given by 

\i,p i J 

( 5 ) 

where 

^ |e,) (e*| -f \gp) {gp\+ujd\d) {d\. (6) 

i p 

Here \gp) are the ground states, \ei) are the excited states, 
and \d) is the distant state. 

The time-dependence of the density matrix p is com¬ 
puted by numerically solving the Liouville equations with 
dispersion in the rotating wave approximation. In partic¬ 


ular, the density matrix and Hamiltonian are expressed 
in terms of the quantities 


PeiQp — Pei 
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where uj^g is the carrier frequency of the laser field cou¬ 
pling the excited and ground states and is the carrier 
frequency of the laser field coupling the excited and dis¬ 
tant states. 

We now make the rotating wave approximation by ne¬ 
glecting all oscillatory terms in the equations of motion 
with frequencies given by a sum of optical frequencies. 
We also include total decay rates T^ and specific channel 
decay rates y^p. The resulting full set of equations of mo¬ 
tion of the density matrix of an Mg + Me system with a 
distant state \d) is given by, 
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Defining N = Mg + Me and taking into account the 
Hermiticity and preserved trace of the density matrix, 
this is a system of equations in N real and {N"^ -I- N)/2 
complex independent variables, as any one of the N + 1 
diagonal density matrix elements can be described as a 
function of only the other N states. Numerical solutions 
to this set of equations are computed using the built-in 
numerical differential equation solver of Wolfram Mathe- 
matica [l^. For a given set of fixed parameters, the com¬ 
putation of the time-dependence of the density matrix for 
a sixteen-level system requires about thirty seconds on a 


moderately fast personal computer (Intel 17-3770 proces¬ 
sor at 3.4 GHz) for a time range of 100/F, where F is the 
smallest of the decay rates F^. 

The behavior of an atom or a molecule with velocity v 
in the lab frame is modeled by making the replacement 
z ^ vt in the spatial dependence of each Rabi frequency. 
This assumes that only motion along the laser beam axis 
is relevant; effects of transverse velocity across an inho- 
mogenous beam profile are not considered here. As in 
previous work [Til m, in the force on the system at 
a particular velocity is then obtained by application of 
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FIG. 2. The 2+1 system simulated in Section ITITl Decay rates 
and ground-state splitting are indicated; the ratio of Rabi fre¬ 
quencies on each transition was set to the ratio of the square 
roots of the branching fractions. 


Ehrenfest’s theorem which gives, for the Hamiltonian in 

Eq.m 


FIG. 3. (Golor online) BCF force from simulations of a 
lamdba-type system, with the BGF carrier frequency reso¬ 
nant with the |gri) —>■ |ei) transition (top) or the \g 2 ) —>■ |ei) 
transition (bottom). Force contours as a function of the total 
Rabi frequency 11*°* and velocity (left) reveal the optimum 
n*°* for each transition (dashed line) which, though at dif¬ 
ferent irradiances, result in similar force vs. velocity profiles 
(right). 


F = -h 


Re 


+VRe[pe.dVH: 


(9) 

Here the VHij are obtained symbolically from Eqs. (THU 
prior to the z ^ vt replacement discussed above. This 
expression gives an immediate general idea of how to 
achieve strong forces: strong coherences must be main¬ 
tained and must change in phase with the driving field. 


III. 2+1 ‘A’ SYSTEMS 


A simple 2+1 lambda-type system without a distant 
dark state was simulated to facilitate systematic study 
of the effects of multiple levels on the bichromatic force. 
The system is set up such that each of the two lower 
states can have a different electric dipole coupling to the 
excited state and a different zero-field energy (Fig. [2]). 

Since each transition has a different coupling to the 
electric field, the BCF laser irradiance produces two dif¬ 
ferent Rabi amplitudes, f^eigi assume in 

this section that each transition is optimally driven and 
that both laser fields have the same polarization, so that 
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( 10 ) 


To aid comparisons, we define a total Rabi frequency 
amplitude H*°* as the quadrature sum of the two Rabi 
amplitudes. 


There are two regimes in which the force resulting 
from BCF illumination is relatively easy to explain, the 
small-detuning regime {S << — Wg^l) and the large- 

detuning regime {6 >> \ujg^ — ujg^\). 

In the small-detuning regime, our simulations show 
that a central BCF frequency resonant with one tran¬ 
sition or the other is preferable to one tuned near the 
weighted average energy of the lower states. When the 
carrier frequency w is on resonance with one of the tran¬ 
sitions, the other transition sees the light as essentially 
a single-frequency laser. Figure [3] shows calculated force 
profiles for resonance with each of the ground-state levels 
|gi) and \g 2 )- While it is necessary to supply more laser 
irradiance on the weaker transition to achieve the optimal 
force at a given detuning, this force is comparable to the 
optimal force attainable on the stronger transition. This 
might be expected at first glance, given that the two-level 
BCF at optimal irradiance depends only on the bichro¬ 
matic detuning. What may confound such an argument 
is the possibility of differing overall coupling strengths 
when the system reaches a periodic quasi-equilibrium. 
Thus we discuss this equilibrium in some detail. 

In a two-level system with resonant monochromatic il¬ 
lumination and without radiative decay, the population 
cycles between being fully in the ground state and fully 
in the excited state. When radiative decay is included, 
the system eventually damps to an equal mixture of the 
ground and excited states. With the addition of third 
state coupled to the excited state, assuming that the 
original transition is resonantly driven, the steady-state 
excited population becomes 


= \/(+.„)" + (11) 


TA^ 
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( 12 ) 
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FIG. 4. Velocity-averaged force for BCF with carrier fre¬ 
quency resonant with the l^i) —>■ |ei) transition of the sys¬ 
tem shown in Fig. O The Rabi frequency amplitude is 
= ■\/3/2 S throughout. Points show results of direct 
simulation and the line shows the two-level force times an es¬ 
timated participating fraction calculated for monochromatic 
excitation (see text). Beyond S « 25r, the system transitions 
to the intermediate-detuning regime. 


where is the detuning of the light from the |ei) O 

[( 72 ) transition. We see from this that for finite Rabi 
frequencies, if the second transition is driven with a 
non-zero detuning, the system remains at least partially 
bright (a full expression can be found in Ref. [^). 
Even though a true coherent dark state does not exist, 
the excited-state population is reduced compared to a 
resonantly-driven two-state system, and some population 
collects in the off-resonant state. For the system shown 
in Fig. [21 if both lower states were coupled to the excited 
state by the same monochromatic light resonant with the 
|ei) O jgi) transition (and therefore = ujg^ — Wg^), 

the combined steady-state population in levels |ei) and 
Igi) would rise with increasing Rabi frequency from zero 
to a maximum of approximately two-thirds that occurs 
at a Rabi frequency beyond the small-detuning regime. 
Likewise for the dichromatic case in the small-detuning 
regime, if the bichromatic carrier frequency is resonant 
with one transition, the other transition is expected to 
experience no appreciable bichromatic force. We expect 
in this case that the bichromatic force in the lamdba-type 
system will be reduced by a factor of the “participating” 
fraction, which can be estimated by adding the popula¬ 
tions in the excited state and resonantly-coupled ground 
state under monochromatic illumination. 

This is verified in Fig. |4l for which the example system 
is simulated using a range of small bichromatic detunings 
6, in each case keeping the carrier uj on resonance with 
the \gi) —>■ |ei) transition and using the optimal Rabi 
frequency = y^3/2(5. In contrast to the two-level 

BCF, the average force is nonlinear in 5, roughly scaling 
as 6^. Up to a bichromatic detuning of one-quarter of the 



Log,g(Bichromatic detuning / ground state splitting) 

FIG. 5. (Color online) Velocity-averaged force in an asymmet¬ 
ric lambda-type system, scaled by the bichromatic detuning 
5, for a range of <5. The contours are given in units of Tik/2. 
Results show a “dead zone” in the force near S = — ojg^ \ 

(vertical dashed line), and very different optimal values of 
in the small- and large-detuning regimes (horizontal 
dashed lines). 


ground-state energy splitting, this force closely follows a 
curve defined by the product of the two-level bichromatic 
force and the participating fraction calculated from the 
from the equilibrium populations for monochromatic ex¬ 
citation. An effective monochromatic Rabi frequency was 
included as a ^-independent fitting parameter in calcu¬ 
lating this participating fraction. The best-fit effective 
Rabi frequency was found to be 2.040(19) 

The full-width at half-maximum (FWHM) of the force 
vs. velocity profiles was found to scale linearly with de¬ 
tuning, as it does for the two-level BCF, at detunings 
S < 25r. The best-fit FWHM is 0.706(12) S/k, somewhat 
wider at the half-maximum point than for two-level sim¬ 
ulations [ 1 ^ . due primarily to a softening of the sharp 
velocity dependence near the edges of the force profile. 

In the intermediate-detuning regime (6 ~ jwg^ — ujg^ |), 
the force drops off significantly (Fig. |5|). Even adjust¬ 
ing the applied laser irradiance does not fully restore the 
force. In this regime, it is impossible to treat the second 
state either as an incoherent perturber or as a coherent 
contributor to the BCF. Clearly this region of the pa¬ 
rameter space should be avoided. 

By contrast, as the detuning increases further and en¬ 
ters the large-detuning regime, the force re-emerges with 
an optimal lower than the optimal value for either 
transition individually. This corresponds to the max¬ 
imum visible in the lower right-hand corner of Fig. [5] 
In this high-detuning regime where 6 >> \ujg^ — 0 Jg 2 \, 
the splitting between the lower states becomes negligi¬ 
ble and a transition is driven between the upper state 
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FIG. 6. (Color online) Force due to large-detuning BCF in the 
Fig.[2]system. Data was smoothed over a velocity range of 0.1 
5/fc. Both the force near zero velocity and the velocity range 
of the force increase with increased bichromatic detuning. 

and a superposition of the two lower states. The to¬ 
tal Rabi frequency then approaches its two-level value of 
= ^/^ 6 . 

Once the dead zone of low force in the intermediate¬ 
detuning regime is exited, the force and velocity range 
both grow linearly with 5, continuing up to quite large 
bichromatic detunings as shown in Fig. [51 The slopes 
for the force and the velocity range are found by least- 
squares fitting to be 0.258(22) Tik/2 and 0.123(21) 1/fc, 
respectively. 

From these results, we can give a general prescription 
for large detunings that the optimal irradiance when mul¬ 
tiple lower states are coupled to a single upper state is 
the one that p rodu ces a total quadrature-summed Rabi 
frequency of y'3/2 S for the transitions “surrounded” by 
the two bichromatic frequencies. In the small-detuning 
regime, where only one transition is addressed, the irra¬ 
diance should instead be chosen s o th at Rabi frequency 
for this individual transition is •\/3/2 6. In both cases, 
we can also say that the force and its velocity range in¬ 
crease monotonically with d. Finally, we have learned 
that the intermediate-detuning region generally leads to 
very small forces, and is to be avoided. 


IV. REPUMPING 

For molecular applications of the BCF, it will often 
be the case that radiative decay losses to distant dark 
states are significant over the desired interaction time, 
requiring a repumping scheme to recover the population 
from these states. To understand the effects of incoher¬ 
ent population loss and subsequent repumping in isola¬ 
tion from other complications, we have simulated a 1-1-1 
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FIG. 7. In a U-l system with 5% of decays going to a distant 
state, the achievable force depends on the Rabi frequency for 
repumping out of the distant state. The force exhibits a clear 
optimum in the repumping Rabi frequency that depends on 
the bichromatic detuning 5. An optimum BGF Rabi ampli¬ 
tude of sj'ij'i. 6 is assumed throughout. Solid lines are smooth 
curves to guide the eye. 

system that also includes a distant dark state. A contin¬ 
uous repumping laser field is used to transfer population 
from the dark state to the same excited state used by the 
cycling BCF transition. Because radiative decay is inco¬ 
herent, the effects of multiple dark states are expected to 
be qualitatively similar to those simulated here for a sin¬ 
gle dark state, apart from the obvious impact of increased 
level degeneracies. 

Figure [7] shows the calculated bichromatic force as a 
function of the repumping laser Rabi frequency for an 
example in which which 5% of decays are to the distant 
state. For each of the several detunings 5 shown in the 
figure, a broad but obvious maximum in the force is ob¬ 
served, indicating that there is an optimum value for the 
repumping laser irradiance. This contrasts starkly with 
more familiar case of repumping for the incoherent radia¬ 
tive force, for which a simple saturation behavior is typi¬ 
cally observed and the repumping irradiance is otherwise 
noncritical. The difference can be attributed primarily 
to decoherence induced when the repumping rate is too 
rapid. 

For the example in Fig. [3 the optimum repumping 
Rabi frequency was found to very closely fit a set of 
curves proportional to and this scaling is expected 

to be generally applicable. Using a least-squares fit, the 
optimal frequency was determined to be 

= 0.536(14)VS^. (13) 

Comparing the BCF with repumping at optimal pa¬ 
rameters to the two-level BCF, the velocity range remains 
unchanged, but the peak force is reduced. For the case 
of 5% decay to a distant state, the force is reduced to 
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FIG. 8. At a fixed bichromatic detuning of lOOF, the op¬ 
timal repump Rabi frequency (top) and the optimized force 
on molecnles with near-zero velocities (bottom) depend on 
the branching ratio of excited state decays to a distant dark 
state. For very small decay rates to the dark state, the force 
approaches a value given by the two-level force multiplied by 
a statistical factor equal to the participating state fraction 
(dashed line, bottom). 


approximately 50% of the two-level force. Further simu¬ 
lations were carried out to compare the reduction in force 
over a range of decay branching ratios, with the bichro¬ 
matic detuning and overall excited-state decay rate re¬ 
maining fixed. As seen in Fig.|5l as the branching fraction 
of spontaneous decays to the distant state increases, the 
peak force decreases monotonically. For very small dark- 
state branching fractions, the force converges to 2/3 of 
the two-level bichromatic force, the value expected based 
on steady-state population statistics. At large branching 
fractions, the peak force decreases approximately linearly 
to zero. The optimal repumping Rabi frequency increases 
with the decay fraction to the dark state over most of the 
range, but begins to decline when more than 70% of the 
decay is to the dark state. 

It should be kept in mind that one of the key advan¬ 
tages of the bichromatic force is the rapid stimulated 
transition rate, which allows a much larger total momen¬ 
tum transfer than incoherent forces if the available inter¬ 
action time is limited. Thus for weak out-of-system decay 
channels, it will often become possible to omit a repump¬ 


ing scheme and allow the loss, as the time it would take 
to accumulate a large population in the distant state may 
be long compared to the needed interaction time. Nev¬ 
ertheless, for the case of rapid decay to dark states the 
need for repumping will be unavoidable. 


V. DARK STATE DESTABILIZATION IN 
MULTILEVEL SYSTEMS 

Here we consider the effects of dark states within the 
transition manifold, which differ fundamentally from the 
distant incoherently-coupled dark states treated in Sec. 
IlYl These dark states can form either when the electric 
dipole selection rules make certain levels inaccessible, or, 
for the case of degenerate levels, when a coherent dark 
state forms in which the excitation dipole amplitude to 
a superposition state coherently cancels. For any sys¬ 
tem in which the ground-state multiplicity of degenerate 
projection quantum states m exceeds the excited-state 
multiplicity by at least two, there will be two dark states 
regardless of t he p olarization of light used to couple the 
two manifolds [^. A number of schemes have been de¬ 
veloped to prevent or reverse optical pumping of the pop¬ 
ulation into these dark states [1, [2l[ , which would oth¬ 
erwise cause the BCF to rapidly diminish to zero. We 
investigate the two most obvious schemes, one involving 
the application of a skewed dc magnetic field and the 
other, rapid switching of the optical polarization state. 


A. Skewed Magnetic Field 

The first method is to apply a dc magnetic field of 
magnitude H at an angle 0be relative to the principal 
quantum axis defined by the optical polarization. For a 
state with angular momentum J and degeneracy 2 J -|- 1, 
this causes a remixing of the mj levels and leads to Zee- 
man shifts proportional (to first order) to mj. This ap¬ 
proach has been analyzed in detail for monochromatic 
excitation in Ref. [2lj and has recently been used suc¬ 
cessfully for radiative cooling and trapping of SrF by the 
DeMille group Here we define the magnetic field to 
lie in the x-z plane, where it can be written in spherical 
tensor notation as 


B = sin(0_B£;)Tj ^+Bcos{9BE)Ti—-^s\n{6BE)Ti^ 

(14) 

The magnetic-field Hamiltonian is modeled by adding 
the additional terms 

^ = MS X! (5pl B ■ {gsS + qeL) \gq) 

+ ^^BY^{e^\B ■ {gsS + gLL)\ej) (15) 
'i'J 
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FIG. 9. (Color online) Contours of velocity-averaged force 
for a ^Po system driven by a linearly-polarized BCF 

field of 5 = lOOr and H = 122.SF, with a dc magnetic field 
of varying amplitude and angle. The dashed line that follows 
the ridge of the force defined by %/2B sin(0) = 12r. 


to the Hamiltonian given in Eq. [S] Magnetic interactions 
with nuclear spins are neglected because they are gener¬ 
ally small, and magnetic coupling between ground and 
excited states is neglected because for small dc magnetic 
fields this coupling is extremely small compared with typ¬ 
ical optical energy level separations. 

In Ref. [2l| the steady-state excited fraction is ex¬ 
amined in a ^Si -O' system driven by linearly- 

polarized monochromatic light with a skewed dc mag¬ 
netic field. We have simulated this same system, replac¬ 
ing the monochromatic optical field with a BCF field. 
We scale the magnetic field units such that the first-order 
Zeeman shift is equal to mjB cos{9be), and express B 
in units of the total upper-state decay rate F. 

For a fixed bichromatic detuning S = lOOF and Rabi 
frequency = -\/3/2 S for the single transition being 
directly driven, simulations with varying magnetic field 
strengths and skew angles reveal a clear optimum force in 
the B-9be plane, evident in Fig. [S] There is also a ridge 
in the force that corresponds to a fixed remixing rate, up 
to a skew angle of about 60°. By fixing B, and Obe at 
their optimum values as determined by these simulations 
and varying ilP, we confirm that a Rabi frequency of 
= -^3/2 S remains optimal in this system. 

At varying detunings, always keeping at its two- 
level optimum, the magnetic field parameter space was 
searched to find how the optimum parameters change 
with detuning. It was found that the optimal intra¬ 
ground-state coupling scales with the square root of the 
detuning and that the optimal level splitting scales lin¬ 
early with the detuning, and we find using least-squares 



1—'—I—'—r 

80 90 100 


Bichromatic detuning (units of F) 


FIG. 10. At optimum magnetic field parameters and Rabi 
frequency, the bichromatic force with magnetic remixing in 
a <-4 ^Po system results in an average force that scales 
linearly with bichromatic detuning and is approximately one- 
third of the ideal two-level bichromatic force. 


fitting that these optimal values are 

V2Bsin{0BE) = 1.184(31)V^, 

Bcosi0BE) = 0.2054(35)5. (16) 

The y/^B sm{0 be) parameter plays the same role here as 
the repumping Rabi frequency in Sec. IlYl which is to re¬ 
introduce lost population to the cycling transition. That 
both have optimum values that scale in the same way as 
•\/^ stands to reason. Eq. (THlcan be rearranged to give 
expressions for the optimum magnetic field strength and 
angle in this system, directly: 


B = 5\/0.700(37)r/5 -h 0.0422(14) 

0be = tan-i (4.07(13) v^) (17) 

If the magnetic field parameters are kept optimized, 
the velocity-averaged bichromatic force remains linear 
with the bichromatic detuning as shown in Fig. [TO] A 
least-squares fits gives the constant of proportionality as 
F/S = 0.1203{27)hk . This is just slightly greater than 
1/3 of the ideal two-level proportionality oi F/5 = hk/ir 
(l/37r = 0.1061). Since the remixing process blends one 
participating ground state with two non-participating 
states, a multiplier close to 1/3 is indeed expected. 

For other multilevel systems the details will differ, but 
the basics of this remixing scheme can be expected to 
be similar, as will the scaling rules. An example of level 
mixing in a much more complicated system is given in 
Section IVl Bl where a detailed 16-level simulation is car¬ 
ried out for the bichromatic force on a molecular beam 
of CaF. 
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B. Polarization Switching 

An alternative method for destabilizing dark states 
is to dynamically switch the polarization of the driv¬ 
ing fields, thereby dynamically switching which state 
is dark. This can generally be accomplished using the 
counterpropagating-beam geometry of a BCF configura¬ 
tion, with the notable exception of four-level transitions 
for which J" = 1 in the lower state and J' = 0 in the 
upper state. While the simple structure of these transi¬ 
tion is often advantageous!^, dark states can be avoided 
only if three polarization states are available, requiring a 
new laser beam at a skewed angle[2lj. This awkward ar¬ 
rangement strongly disfavors the polarization-switching 
method relative to the alternative of using a skewed B 
field. 

We consider transitions of the type J' = J" + 1 with 
J" > 1. Here the dark states can be avoided by alter¬ 
nating between and cr~ polarization. For example, 
in a J" = 2 J' = 1 transition, defining the quanti¬ 
zation axis along the fc-vector of the stimulating light, 
cr'*' light leaves m" = 1 and m" = 2 dark, while cr“ 
light leaves m" = — 1 and m" = —2 dark. Clearly the 
dark-state spaces for these polarizations do not overlap, 
and the necessary polarization switching can readily be 
performed on a collinear pair of beams. 

To test whether polarization switching during BCF 
still allows significant forces, a ^82 ^ ^Pi system was 
simulated with the polarization of the BCF field switched 
between cr''" and a~. Since the pulse area of the bichro- 
matic beat notes is a critical BCF parameter, a sinusoidal 
variation between polarization states can be expected to 
give poor results because it leads to non-constant pulse 
areas. Instead, we simulate a simple on-off modulation 
pattern in which each polarization is kept on for a time 
T before switching to the other polarization. 

During each step in this cycle, population will gradu¬ 
ally be optically pumped into the two dark mj states 
that cannot be re-excited, diminishing the force com- 
mensurately. Each switch in polarization constitutes a 
rephasing of the system after which coherence on the new 
transition (and thus the optical force) will gradually be¬ 
gin building up. An optimal switching period is therefore 
expected that balances these two effects. 

The initial phase of the bichromatic optical field influ¬ 
ences how quickly an incoherent population equilibrates 
under bichromatic driving. In most of our simulations, 
this phase is set to 9 = 45°, which gives the quick¬ 
est equilibration, and the transient behavior is allowed 
to die away before evaluating the force. However, this 
transient behavior takes up a significant portion of each 
polarization-switching cycle and must now be treated 
more systematically. Simulations for which the global 
phase is reset to a particular value 9 at each switch of 
polarization show that 9 = 45° does indeed result in the 
strongest force, and 0 = 0° the weakest, actually exhibit¬ 
ing a weak force in the reverse direction. 

The ^82 O ^Pi system has a further complication be- 



FIG. 11. Bichromatic force for a ^82 -H- transition as 
a function of the switching period T between and g~ 
polarization states. The Rabi frequency is set to its opti¬ 
mum value for the \m"\ = 1 m' = 0 transition or for the 
\m"\ = 2 |m'| = 1 transition, giving the two curves above, 

and the bichromatic detuning is lOOF throughout. 


cause the three transitions driven by cr-polarized light are 
not equally strong. Therefore it is not immediately ap¬ 
parent what laser irradiance is optimal. Simulating this 
polarization switching scheme at 5 = lOOF over a range 
of switching times T showed that a maximum in the force 
is achieved when setting the Rabi frequency to the value 
optimal for the |m"| = 1 o m' = 0 transition. Also, 
a switching time of T ft! 5/r is found to be optimal as 
shown in Fig. [TTJ with a shallower fall-off on the long-T 
side (Fig.) than for short T.. The peak force in these 
simulations is approximately 207ifcF/2, which is about 
1/3 of the peak two-level force at the same bichromatic 
detuning, 64?iA:r/2. 

If the switching cycle is not phase-locked to the global 
bichromatic phase, the random starting phase at each po¬ 
larization switch will cause a reduction in the force com¬ 
pared with these simulations. As expected,the polariza¬ 
tion switching period must be at least comparable to the 
upper state decay rate (^ F/5 in these simulations) and 
the transition time between polarization states should be 
short compared to l/i5, typically requiring nanosecond- 
scale switching. 

For comparison, dark-state destabilization by a skewed 
magnetic field was also simulated for the same ^82 
^Pi model system, with 7r-polarized BCF helds. With 
the Rabi frequency optimized for |m = 1| transitions, 
such that two of the three available transitions were being 
driven with optimal Rabi frequency, the peak force over a 
considerable range of magnetic field parameters is found 
to be 24hkr /2. Unlike the polarization switching scheme, 
no fine-tuning of phases or timing is required. At least in 
this system, magnetic field remixing is the better choice, 
but if a particular application cannot tolerate external 
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mp=-2 mp=-1 mp=-0 mp=1 mp=2 



showing zero-field energy spacings. The hyperfine splitting of 
the B levels is taken to be negligible compared to the X state. 
The allowed transitions for yr-polarized light are indicated. 

magnetic fields, polarization switching remains a viable 
alternative. 


TABLE I. Electric dipole matrix elements for the B X 
Pii(0.5)/^(3i 2(1.5) branch in CaF, in units of the total B o 
X dipole transition moment. 




o 

II 

fe, 


R'=l 


F” 

/f 

rrip 

m'p^O 

f 1 

mp=—l 

O 

II 

s 

/ 1 
mF = l 


-2 

0 

-0.5774 

0 

0 


-1 

0 

0.4082 

-0.4082 

0 

2 

0 

0 

-0.2357 

0.4714 

-0.2357 


1 

0 

0 

-0.4082 

0.4082 


2 

0 

0 

0 

-0.5774 


-1 

-0.5743 

-0.0421 

-0.0421 

0 

1+ 

0 

0.5743 

0.0421 

0 

-0.0421 


1 

-0.5743 

0 

0.0421 

0.0421 


-1 

0.0595 

-0.4061 

-0.4061 

0 

1" 

0 

-0.0595 

0.4061 

0 

-0.4061 


1 

0.0595 

0 

0.4061 

0.4061 

0 

0 

0 

0.3333 

0.3333 

0.3333 


VI. CALCIUM MONOFLUORIDE 
A. Structure and Parameters 

With an eye towards experimental realization of the 
BCF using calcium monofluoride, the structure of the 
(0 — 0 ) Pii(1 .5)/'^(5 i2 (0.5) branch was 
examined. This scheme differs in some details from the 
A X transition that we have previously discussed 
semi-quantitatively in Ref. although the differences 
are for the most part minor. The fine and hyperfine sub- 
levels for this transition are shown in Fig. [T2l In the X 
state 7V=1 manifold, J is a poor quantum number due 
to mixing of the two F=N=1 levels. The two eigenfunc¬ 
tions \F = 1=*=) can be written in terms of the J=\j2 and 
J=3/2 levels with F=1 by introducing a mixing angle 0, 

/|F=1+,mi.)\ ^ / cos(/> sin(/i\ /|3/2,mi.)\ 

V|T=l“,mF)/ y—sin^ cos4>) \\l/2,mF)] 


We define a relative dipole transition element between 
states I*) and \j) as 


E,{ATy\d)\j) 

(A', 5, S' = 1/2| I A, 5, S = 1/2) ■ 


( 20 ) 


We calculate these relative dipole transition matrix ele¬ 
ments, following Ref. [1^ with the exception that A A = 0 
rather than ±1, and then introduce the J-mixing dis¬ 
cussed above to obtain the relative dipole transition ma¬ 
trix elements between the eigenfunctions (Table |I|. 

The normalization factor can be estimated by working 
back from the measured radiative decay lifetime of the B 
state, TB = 25.1 ± 4 ns [13. The characteristic lifetime 
of a state before it decays along a particular transition is 
given by the inverse of the Einstein A coefficient, 


Bj — 1/^ij 


ShX^So 
167r3 \dij\'^ 


( 21 ) 


The mixing angle was determined to be 41.18° via the 
methods of Ref. [l^, in which the structurally similar 
YbF molecule was examined. We used the molecular 
constants for CaF listed in Ref. [^ . 

When the pure J states are expressed in case 
{ap) notation, | A, S', S, O, J,/, F, Mi?), the electric 
dipole transition matrix elements can be calcu¬ 
lated as in Ref. [13 up to a shared factor 
(A', S, E' = 1/2| T^^\d) |A, S, S = 1/2). Additionally, 
the squared sum of matrix elements is normalized to this 
factor. 


E E 

P J",F",Mf 


(J',F',M/,|r«(d)|J",F",M") 


(A',S,l/2|To(')(d)|A,S,l/2) 


(19) 


The total lifetime of a state is the inverse of the sum 
of the A coefficients for all possible decays. Outside of 
the channels we are considering, other rotational decays 
on R o A are forbidden, but other vibrational decays 
are allowed, as are decays to other electronic states. The 
Franck-Condon factor for the R o A (0-0) band has 
been calculated to be 0.999 [3l|. The decay branching 
ratios to different electronic states scale as |<iy The 
dipole transition moments have been estimated via lig¬ 
and field modeling [13 and the transition frequencies are 
known [ 13 , so the ratio of decay rates to the can be es¬ 
timated: 

= (1.73 eao)^(18833cm-i)3 : (0.58 eao)^(2307cm“^)3 
= 4840 : 1. (22) 
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TABLE II. g-factors for the six F levels involved in the B -f-s- 
X Pii(0.5)/^Qi 2(1.5) branch in CaF, scaled by the electron 


spin g- 

Factor gs. 



F''=Q F''=l- 

F"=i+ F"^2 F'^0 F'=l 

gp/g!, 

0.0 0.320 

-0.070 0.25 0.0 0.5 


Thus, approximately 99.98% of radiative decays from the 
B state produce molecules in the X state. Multiplying 
the electronic, vibrational and rotational factors, 99.8% 
of decays from B, u=0, A^=0 produce the X, u=0, iV=l 
state. 

This allows an estimate of the total A coefficient for 
the B X Pii( 0 . 5 )/^( 5 i 2 ( 1 . 5 ) branch. 


A = 


0.998 

tb 


16 tt^ 


(A', S, 1/21 T(^\d)lA, 5,1/2} 


SheoX 


(23) 

where A = 530.9 ±0.1 nm is the transition wavelength. 
Rearranging, we can estimate that the matrix element is 


(A',5,l/2|TW(d)|A,5,l/2) =1.7(l)eao, (24) 


which is consistent with the value from Ref. [s^ of 1.73 
eap. Combining Eqs. |4] and [20] with the irradiance of a 
monochromatic plane wave gives practical units for the 
Rabi frequency amplitude of a transition in this chan¬ 
nel for a given per-bichromatic-component irradiance I 
expressed in W/cm^, 


= KijV7(27r X 60 MHz), (25) 

assuming that the light is optimally polarized to drive 
the transition. 

Assuming that magnetic remixing of dark states is 
used, the response to a dc magnetic field must be un¬ 
derstood as well. The magnetic Hamiltonian in Eq. [121 
can be evaluated using equations (8.183) and (8.184) of 
Ref. [S^, with a few obvious modifications to allow for 
the spherical components p 0. These equations give 
the complete picture, but the six F states can be gen¬ 
erally characterized in a more intuitive way via their g- 
factors gp, which for the two F" = 1 levels must take the 
J-mixing of the eigenstates into account (Table [III). 


B. Simulations 

Previous treatments of the bichromatic force in mul¬ 
tilevel systems i,[i3 have been restricted to analyses 
based on two-level calculations combined with statistical 
arguments. These treatments cannot account for multi¬ 
level coherent effects such as those seen in Sec. m which 
were included along with detuning and line-strength con¬ 
siderations in a poorly-determined “force reduction fac¬ 
tor,” T]. Examining the structure of the ^11“^ 


F' = 2, +48.8 MHz 
F" = r, +23.9 MHz 

F" = 0, -22.9 MHz 


F" = r. -97.6 MHz 


F' = 0 F' = 1 F' = 1 F' = 1 

mp = 0 /Hp = -1 mp = 0 mp=1 



FIG. 13. The four subsystems defined by isolating the Amp = 
0 transitions of the B -(P)- X Pii(1.5)/^(5i2(0.5) branches in 
CaF. Relative dipole transition matrix elements are indicated 
for each transition. 


(0 — 0) Pii(1 .5 )/^Qi2 (0.5) transition in CaF (Fig. [T2|), 
it is easy to see that it cannot be fully decomposed into 
two-level systems even for the simplest case of 7r-polarized 
light. Each excited-state sublevel is coupled to either two 
or three ground-state levels, and thus the coherent effect 
contribution to 77 cannot be neglected. 

However, Fig. [T^] also suggests that some simplifica¬ 
tion may be possible without giving up much accuracy. 
No ground-state sublevel is laser-coupled to more than 
one excited-state sublevel, so that the entire system can 
be viewed as two 2±1 and two 3±1 subsystems which 
are weakly coupled to each other by incoherent radia¬ 
tive decays, in addition to two dark ground-state levels. 
We start by taking this simplified approach and model¬ 
ing these four subsystems independently, neglecting the 
radiative coupling between them. Because the coherence 
terms of the density matrix decay at a rate set by the to¬ 
tal excited-state decay rates, we proportionally increase 
the decay rates within each subsystem to compensate for 
the neglected decays. 

The levels in the X, N" = 1 states span 146.4 MHz, or 
23.1 r in units of the excited-state decay rate. As seen in 
Sec.[TTTl so long as the bichromatic detuning is larger than 
the energy splitting, the optimal 17*°* will be ^J‘^l/2 5. 
Within each of the four subsystems, the quadrature sum 
of the Kij is l/'x/S, so the same laser irradiance can drive 
every one of the subsystems with the same 17*°*. 

Using the optimal 17*°* at each of <5 = SOT, 50r and 
lOOr, each of the four subsystems was independently 
simulated for a range of BCE carrier frequencies that 
spanned the entire set of transitions. The results from 
each subsystem were then combined together with a 
weighting given by the fraction of ground state sublevels 
included in each subsystem. This weighting is based on 
the assumption that the system begins with a statistical 
distribution of population in the ground states and that 
each subsystem will maintain that population through¬ 
out the interaction with the BCE fields. As shown in 
Fig. 1141 the averaged force on molecules with near-zero 
velocities was found to peak at a carrier frequency de- 

































12 


12 - 


. 10 


£ 8 


> 

CD 4 

'o 

o 2 

a> 

> 






5=30r 

■<5=50r 

<y= loor 


A-^'A' 








-150 -100 -50 0 50 100 150 

Carrier frequency detuning from com (MHz) 



FIG. 14. A weighted average of the force on isolated subsys¬ 
tems of the B X Pii(1.5)/^(5i 2(0.5) branches in CaF gives 
an estimate for the optimal carrier frequency detuning from 
the transition center of mass of the system and for the peak 
force. Resonances are indicated with vertical dashed lines. 


FIG. 15. Simulations of the dichromatic force for the sixteen- 
level B X Pii(0.5 )/^Qi 2(1.5) system in CaF, with dark- 
state remixing by a dc magnetic field. These forces far ex¬ 
ceed the radiative force of ~0.024 hkV/2 realized in Ref. [^. 
The strength and velocity range of the force increase with the 
dichromatic detuning in the range of detunings simulated. 


tuned somewhat from the hyperfine center of mass, at a 
location between the F" = 2 and F" = 1+ components. 
Transitions from these two levels were the strongest tran¬ 
sitions in each of the four subsystems. 

Using optimal parameter values, estimated values for 
the key properties relevant to longitudinal slowing of 
a molecular beam are summarized in Table Ell in the 
columns labeled as Method (b). The table also includes 
estimates of the these same quantities based on the semi- 
quantitative statistical method of Ref. @ , where applica¬ 
ble, as well as estimates based on the complete 16-level 
simulation that we describe next. The subsystem ap¬ 
proach must be regarded as an approximation both be¬ 
cause it relies on assumptions about the population bal¬ 
ance, and because it cannot include the effects of remix¬ 
ing dark states. Instead it is assumed that on aggregate, 
a good remixing scheme will keep the dark populations 
at about the statistical fraction determined by the total 
number of states. The results of Sec. |V] support this esti¬ 
mate, but the details are revealed only by modeling the 
full system. 

Using the optimal laser irradiance determined from 
subsystem-based modeling at each dichromatic detuning, 
the entire set of sixteen levels was then simulated for tt- 
polarized BCF fields and a skewed dc magnetic mixing 
field, using a two-step process. First, the remaining pa¬ 
rameter space consisting of the magnitude and angle of 
the magnetic field and the carrier frequency detuning was 
surveyed at each of several small molecular velocities, in 
order to find optimal parameters. Then, using those pa¬ 
rameters, simulations across the full velocity range were 
carried out. The resultant force profiles are shown in 
Fig. [12] for each of three representative dichromatic de¬ 
tunings, and the properties relevant to longitudinal slow¬ 
ing of a molecular beam are listed in the columns of Ta¬ 


ble Hill labeled as Method (c). 

As seen previously for 2-1-1 systems in Sec. Illll for de¬ 
tunings in the lower portion of the large-detuning regime, 
the force increases slightly faster than linearly with S. 
Comparing to the results from the subsystem approach, 
it is clear that the subsystem approach consistently over¬ 
estimates the peak force and the average excited state 
fraction, and it underestimates the velocity range of the 
force. However, all of these errors are small, with less 
than 10% error in the force and less than 25% error in 
the velocity range across all three simulated detunings. 
In this example the subsystem-based calculations are not 
only simpler, but they also required only 1/5 of the com¬ 
putational time needed for the full sixteen-level calcula¬ 
tion. These considerations may make the simpler and 
more approximate subsystem method preferable for cal¬ 
culations intended for planning purposes. 

Compared to both of our multilevel models, the statis¬ 
tical approach of Ref. @ is, of course, much less accurate. 
It slightly underestimates the velocity range of the force 
at 5 = 30r but this crosses over to a large overestimate 
at (5 = lOOF. This can be attributed to the approxima¬ 
tion that Av = 5/k that does not accont for destructive 
coherent effects that limit the range of the force in the 
intermediate detuning regime. There is also a substantial 
overall uncertainty due to the “force reduction factor” rj, 
which essentially lumps together all of the effects not in¬ 
cluded in the effective two-level model. Rough estimates 
in the spirit of Ref. @ predict that rj should be about 
0.57 at (5 = 30r, increasing to about 0.8 at lOOF. By 
comparing the estimated force to our full 16-level simu¬ 
lations, we can empirically evaluate the actual value of 
r] in each case, giving r] = 0.16, 0.21, and 0.26 for 6 = 
30r, 50F, and lOOF respectively. This verifies that the 
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TABLE III. A summary of results for three different methods of estimating the BCF for vr-polarized light in the sixteen- 
level B X Pii(0.5)/^Qi 2(1.5) system in CaF, assuming no repumping. The methods are (a) the statistics-based estimate 
described in Ref. [^, (b) a weighted subsystem-based multilevel simulation, and (c) a full 16-level simulation of the complete 
system, including dark-state remixing by a dc magnetic field. Some results for Method (a) includes a “force reduction factor” 
rj as in Ref. [^, discussed further in the text. In the final row of the table, we consider the fraction of molecules remaining after 
slowing by 60 m/s, taking into account decays into distant dark states. 



(a) 

^=3or 

(b) 

(c) 

(a) 

5=50r 

(b) 

(c) 

(a) 

5=ioor 

(b) 

(c) 

Irradiance (W/cm^) 

45.2 

45.2 

45.2 

125.6 

125.6 

125.6 

503 

503 

503 

Carrier detuning from c.o.m. (MHz) 

N/A 

-32 

-38 

N/A 

-30 

-38 

N/A 

-47 

-38 

Magnetic field magnitude (Gauss) 

— 

— 

29.2 

— 

— 

29.4 

— 

— 

37.1 

Magnetic field angle 

— 

— 

71° 

— 

— 

61° 

— 

— 

66° 

Force near zero velocity {hkV /2) 

16.477 

2.8 

2.6 

27.377 

6.3 

5.8 

54.677 

15.5 

14.2 

Velocity range (m/s) 

100 

no 

120 

170 

130 

170 

340 

190 

250 

Average excited fraction 

0.14 

0.23 

0.18 

0.14 

0.21 

0.15 

0.14 

0.18 

0.14 

Time to slow by 60 m/s {fis) 

14/77 

84 

91 

8 . 7/77 

38 

41 

4 . 3/77 

15 

17 

Fully-slowed population fraction 

— 

0.21 

0.27 

— 

0.53 

0.61 

— 

0.80 

0.83 


statistical approach is useful for making rapid order-of- 
magnitude estimates, and that as expected the value of 
r] gradually increases with detuning. However, it also 
verifies that the actual numbers are not reliable beyond 
this level of approximation. The method is useful for 
determining feasibility, but is not necessarily reliable for 
determining detailed experimental designs. 

We conclude by comparing our the multilevel BCF 
with the radiative force, working from our sixteen-level 
simulation results. Even for the lowest detuning simu¬ 
lated, the calculated dichromatic force is already much 
greater than the largest radiative force achieved in CaF 
to date. Using the same X-state level in an A o X 
rotationally-closed transition, a radiative force with a 
peak of 0.024(6) fikV/2 has been realized Q. Our re¬ 
sult of 2.6 TikT /2 at d = 30F is larger by more than two 
orders of magnitude. This implicitly assumes that the 
BCF scheme does not require the vio repumping laser 
used in Ref. Because the decay branching fraction 
to dark states is low, the results of Sec. IIVI indicate that 
even with repumping, the force would be reduced only 
by a statistical factor of at most a single order of magni¬ 
tude. However, we now take a closer look to demonstrate 
that the BCF acts rapidly enough to make repumping un¬ 
necessary for a practical decelerator that operates on a 
buffer-gas-cooled molecular beam. 

In Sec. IVI Al we estimate that 0.2% of radiative de¬ 
cays populate non-cycling states. The characteristic out- 
of-system decay time is this tb/{0 -002Pe), where Pe is 
the average excited-state fraction. We find Pe by averag¬ 
ing the calculated excited-state fraction over the velocity 
range in which the force is large. The result is Pe ~ 1/6, 
with variations of just a few percent between the three 
simulations in Fig. [151 This value is also quite close to 
the estimate of 1/7 provided by the statistical approach 
of Ref. 0 . Thus we can estimate an out-of-system decay 
time of about 75 ^s. 

In recent work elsewhere, radiative slowing was used 
to slow a cryogenic beam of SrF by 60 m/s, and this 


proved sufficient to load a molecular MOT [^. A rea¬ 
sonable objective for CaF is to slow it by a comparable 
amount. Using the bichromatic force with a detuning of 
6 = 30F, we estimate that 27% of the illuminated CaF 
molecules with N"=l would be slowed by 60 m/s with¬ 
out repumping, and this percentage increases to 83% for 
6 = lOOF. These results clearly indicate that as long 
as the required laser irradiance is available to drive the 
BCF with detunings of at least 30r, a repump-free op¬ 
tical decelerator is feasible for CaF and other molecules 
with similar near-cycling transitions. 


VII. CONCLUSIONS 

If the coherent optical bichromatic force is to be ap¬ 
plied to molecules or complex atoms, it is essential to 
prove that it can work for an intricately coupled mul¬ 
tilevel system, and to have a reliable method for mod¬ 
eling such systems so that an optimal BCF scheme can 
be devised. We have demonstrated that direct numerical 
integration of the density matrix provides stable and use¬ 
ful solutions for systems with up to 16 sublevels, using 
only the computational capacity of an ordinary personal 
computer. 

For three-level lambda-type systems, we have demon¬ 
strated and quantitatively modeled a nonlinear depen¬ 
dence of the bichromatic force on detuning at small de¬ 
tunings. We have also shown that intra-ground-state co¬ 
herences play a vital role at large detunings, allowing 
lower laser power requirements than would be estimated 
from individual transitions, with the general prescription 
that the quadrature sum of the Rabi frequencies of the 
addressed transitions should be equal to a/ 3/2(5. Fur¬ 
ther, we have shown that the region of parameter space 
where the bichromatic detuning is roughly equal to the 
transition energy splittings is to be avoided, as it results 
in severe reductions in the achievable force. 

When an additional distant dark state is introduced. 
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we show that a repumping laser can be effective, al¬ 
though the requirements for the laser power are more 
specific than for repumping with the ordinary incoher¬ 
ent radiative force. For transitions involving more than 
three sublevels, there are additional complications due 
to the formation of dark states within the system, but 
we show that this can be addressed by adding a small 
skewed magnetic field or by rapidly switching the optical 
polarization state. We have found that the magnetic field 
method is generally preferable. 

We have shown that for many-level systems, a simu¬ 
lation approach in which the system is decomposed into 
individually-evaluated coherent subsystems provides sig¬ 
nificant computational speed-up with associated errors of 
less than 10% in the predicted force. 

Most importantly, we have quantitatively verified the 
predictions of our 2011 paper, clearly demonstrating that 
BCF slowing and cooling of molecules should be both 
practical and promising. Examining the specific example 
of CaF, we show that the required laser power is readily 
attainable and the magnitude of the force is over two 
orders of magnitude larger than the radiative force for 


the same system. Further, the momentum transfer prior 
to out-of-system decay is large enough to allow complete 
deceleration of a cryogenic beam of CaF without the use 
of repumping lasers and over a distance scale of just a 
single centimeter. 

Similar gains relative to the radiative force can be ex¬ 
pected for other small molecules. In general, limited 
momentum transfer sufficient to produce state-selective 
transverse deflection should be feasible for almost any 
small molecule. By contrast, substantial longitudinal 
slowing requires a near-cycling transition because the re¬ 
quired momentum transfer involves thousands of pho¬ 
tons, but the requirements are still much relaxed com¬ 
pared with the conventional radiative force. 
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